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White dwarfs represent the last stage of evolution of stars with mass less than about eight 
times that of the Sun and, like other stars, are often found in binaries!? If the orbital period 
of the binary is short enough, energy losses from gravitational-wave radiation can shrink the 
orbit until the two white dwarfs come into contact and merge. Depending on the component 
masses, the merger can lead to a supernova of type Ia or result in a massive white dwarf! 
In the latter case, the white dwarf remnant is expected to be highly magnetised*!4 because of 
the strong magnetic dynamo that should arise during the merger, and be rapidly spinning 
from the conservation of the orbital angular momentum, Here we report observations of 
a white dwarf, ZTF J190132.9+145808.7, that exhibits these properties, but to an extreme: 
a rotation period of 6.94 minutes, a magnetic field ranging between 600 megagauss and 900 
megagauss over its surface, and a stellar radius of 2, 140* 199 km, slightly larger than the 
radius of the Moon. Such a small radius implies that the star's mass is close to the maximum 
white-dwarf mass, or Chandrasekhar mass. ZTF J190132.9+145808.7 is likely to be cooling 


through the Urca processes (neutrino emission from electron capture on sodium) because of 


the high densities reached in its core. 


Using the The Zwicky Transient Facility! (ZTF) we searched for short period objects that 
lie below the main white dwarf cooling sequence in the Gai?! colour-magnitude diagram (see 
Fig. [1]. ZTF J190132.9+145808.7 (hereafter ZTF J1901+1458) showed promising small photo- 
metric variations. Followup with CHIMERA“, a high-speed imaging photometer on the 200-inch 
Hale telescope, confirmed a period of 6.94 minutes (see Fig. [2). The period of ZTF J1901+1458 


is unusually short for a white dwarf, as white dwarf rotational periods typically are upwards of 


hours! The period and ephemeris are listed in Table|1| 


We undertook phase-resolved spectroscopy using the Low Resolution Imaging Spectrometer 
(LRIS)2 on the 10-m W. M. Keck I Telescope. As can be seen from Fig. 3| the phase-averaged 
spectrum exhibits broad and shallow features that we identify as hydrogen absorption lines in a 
high magnetic field. The presence of a strong magnetic field results in splitting and proportional 
shifting of the zero-field energy levels, leading to line broadening. To identify the field strength, 
we considered all the allowed bound-bound hydrogen transitions (tabulated in) and, as shown in 
Fig.|3, we find a satisfactory explanation for the spectrum. The identified spectral lines are listed 


in Extended Data Table P] 


We find that most of the spectral features are well characterised by a magnetic field strength 
of 800 million Gauss (MG, red horizontal line in Fig. B), comparable to the field detected on 
the most magnetic white dwarfs known?! As the absorption features indicate an average field 
strength over the surface, the field at the magnetic pole is bound to be higher. From the phase- 
resolved spectra (Extended Data Fig. [4] and [5) we see that some of the features become narrower 
or broader depending on the phase, and the feature at ~ 4600 A shifts in wavelength, which 
accounts for the dip at ~ 4500 A in the co-added spectrum. This means that in some regions of 
the surface the magnetic field is as low as 600 MG. The explanation for the photometric variation, 
confirmed by the variations of absorption features with phase, is thus the combination of magnetic 
dichroism and rotation: the high magnetic field causes variations in the continuum opacities and in 


the surface temperature, and therefore, as the star rotates, we detect changes in flux as a function 
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Figure 1: Gaia colour-magnitude diagram. Gaia color-magnitude diagram for the white dwarfs 
that are within 100 pc from Earth and within the SDSS footprint! where the x-axis depicts the 
difference between the Gaia BP and RP bands, and the y-axis the absolute magnitude in the Gaia 
G filter. Solid lines show theoretical cooling tracks for white dwarfs with masses between 0.6 Mo 
(top) and 1.28 Mo (bottom), equally spaced in mass; the atmosphere is assumed to be hydrogen- 
dominated?! and the interior composition to be carbon-oxy geri?! for M « 1.1 Mo and oxygen- 
neon? for M > 1.1 Mọ. Coloured markers indicate white dwarfs for which a magnetic field was 
measured!) ZTF J1901+1458 is shown as a red star, and its location in the colour-magnitude di- 
agram reveals its high mass. Vertical dashed lines indicate the location of the ZZ-Ceti instability 
strip?! (the pre-white dwarf, or DOV, instability strip lies above the plo£283) . Reddening cor- 
rections were applied only to ZTF J1901+1458; as the objects in the sample are close, reddening 
is expected to be small. lo error bars are id than the size of the coloured markers,and are 


omitted for the black background dots for clarity. 
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Figure 2: ZTF J1901+1458 lightcurve The left panels show the binned CHIMERA lightcurve 
phase-folded at a period of 6.94 minutes in the g’-band (a) and in the r'-band (b). The flux has 
been normalised to the minimum of the lightcurve in each band. The amplitude of the photometric 
variation is higher in the g’-band (about 3% peak-to-peak) than in the r’-band (~ 1.596). Addi- 
tionally, the two filters show a difference in phase: the red lags the green by about 51 s. The right 
panels show the similarly normalised ZTF discovery lightcurve in the ZTF g-band (c) and r—band 


(d). The error bars indicate lo errors. 


of the field strength across the stellar surface. Depending on the magnetic field configuration, 
this dichroism can account for up to 10% photometric variation) so it can easily account for the 
3% amplitude observed in ZTF J1901--1458. ZTF J1901+1458’s period could also be consistent 
with non-radial pulsations; however, its temperature and surface gravity place it far away from 
theoretical predictions for known instabilities (see Fig. (I), and its magnetic field may be strong 


enough to suppress gravity-mode pulsations (see the Method section for further discussion). 


To determine the temperature and radius of the white dwarf, we obtained photometric mea- 
surements in the near UV using the UVOT (Ultraviolet/Optical Telescope) instrumenti? on the 
Neil Gehrels Swift Observatory while the Pan-STARRS Survey!’ (PS1) and the Gaia mission 
provided optical photometry and parallax, respectively. We estimated the temperature, reddening 
and radius from the photometry, by comparing the observations with non-magnetic white dwarf 
atmosphere models. Because the white dwarf is very hot, the photometric constraints on the tem- 
perature are weak; however, the precise distance measurement from Gaia allows us to obtain a 
good estimate for the radius. We found the effective temperature, stellar radius and interstellar 
reddening to be Tyg = 46, 000* 2500. K, R = 2, 1404389 km and E(B — V) = 0.044 0015, respec- 
tively (see the Methods section for more details). The radius is smaller than those measured for 
other white dwarfs and only slightly larger than that of the Moon. As explained below, the small 


radius also means that ZTF J1901+1458 maybe the most massive white dwarf yet discovered. 


The mass can be inferred from the mass-radius relation which, as can be seen in Fig. 


is composition dependent. White dwarf descendants of single stars with mass above 1.1 solar 
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Figure 3: ZTF J1901+1458 optical spectrum. The LRIS phase-averaged spectrum of 
ZTF J1901+1458 is shown in black in the upper panels for the blue (a) and red (b) sections. 
In (c)-(d) the spectrum is compared to predicted line positions of Ho, H@ and Hy as a function 
of magnetic field?! showing that the white dwarf is characterised by a field strength of about 800 
MG (red horizontal solid line). The identified absorption features are highlighted by the green 


dash-and-dotted lines and numbered, and the respective transitions are listed in Extended Data 


Table 


masses (Mo) are expected to be mostly oxygen and neon, with traces of carbon, sodium and 
magnesiuni!?P9| Even if born from a merger, compressional heating due to rapid accretion is 
expected to ignite off-center carbon burning?! resulting in an O/Ne white dwarf! We conclude 
that, depending on the composition, ZTFJ1901+1458 has a mass between 1.327 and 1.365 solar 


masses. 


Two other extremely massive white dwarfs for which a radius has been measured are RE J0317- 
8532land WD 18324089231 both with a radius of about 2,500 km. Curiously, they are both variable 
with short periods (6 and 12 minutes, respectively). While WD 18324089 does not show evidence 
for magnetism, RE J0317-853 appears to have a field that varies between 185 and 425 MG over its 
surface! Thus both ZTF J1901+1458 and RE J0317-853, with their rapid rotation, high mass and 


magnetism, are likely remnants of a white dwarf merger. 


At the densities reached in the centre of ZTF J1901+1458, the nuclei of some elements may 
undergo electron capture (also called inverse beta-decay), removing electrons that contribute to the 
degeneracy pressure that keeps the star from collapsing. This lowers the maximum mass that can 
be sustained against gravity and reduces the equilibrium radius for a fixed mass (see Fig. 4| solid 
red curve). If ZTF J1901+1458 has an oxygen-neon internal composition (as is expected from its 
masdi), its central density is right at the threshold for electron capture on ??Na, and its mass is 


within 246 of the highest possible mass for a white dwarf. 


Neutrinos produced by inverse beta-decay escape and carry away energy, contributing to 


the cooling of the white dwarf. Based on the luminosity of the white dwarf?! we estimate the 
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Figure 4: Mass-radius relation. The white dwarf mass-radius relation at the temperature of 
ZTF J1901+1458 (~46,000 K) is shown for different homogeneous compositions neglecting the 
effect of electron captures: carbon C; oxygen O; neon Ne; and the mixture that is thought to result 
from carbon burning (M;), specifically 58% oxygen, 30% neon, 2% carbon, 5% sodium, and 5% 
magnesium?! The lowermost red solid curve (M2) traces the result from the carbon-burning-ash 
mixture including the effect of electron captures (see text). The black marks along each track in- 
dicate the radius at which the central density reaches the threshold for inverse beta-decay of neon 
(circle), magnesium (square) and sodium (diamond). The vertical blue-shaded area indicates the 
observed radius of ZTF J1901+1458 within lo errors. The red triangle indicates the maximum 
mass of 1.354 Mo at a radius of 1,600 km along the track of the mixed composition. The horizon- 
tal black solid lines highlight the mass determined using a radius of 1,910 km for the upper line 


(pure C), 1.365 Mo and of 2,300 km for the low&r line (carbon-burning-ash mixture), 1.327 Mo. 


temperature in the core to be about 2 — 3 x 10" K. At such a high central temperature and density, 
the neutrino cooling of ZTF J1901+1458 will dominated by the "Urca" process?4P3lacting on Na. 
This unusual neutrino cooling makes an age determination difficult. A recent work? studied the 
evolution of Urca-cooling white dwarfs; from the same models, we estimate the cooling age of 


ZTF J1901+1458 to be between 10 and 100 Myr. 


The mass-radius relation in Fig. |4| was calculated assuming that the core composition is 
homogeneous — a good assumption since ZTF J1901+1458 is less than 100 Myr old. However, 
over a few hundred million years, the heaviest elements, including Na, will gradually sink to the 
centré??| If the star lies at the small end of the radius constraint and if at least sixty percent of the 
?3Na manages to sink and undergo beta decay before the core crystallises and sedimentation stops, 
electron-capture on ?"Mg would ensue. The star would shrink and the internal pressure would no 
longer be able to support the star, as the maximum allowed mass for the new composition would be 
lower than the mass of the white dwarf (see the Methods section for a more detailed discussion). 
The star would therefore collapse and heat up, leading to the onset of electron capture onto Ne 
and to the ignition of oxygen nuclear burning. The white dwarf would then undergo a disruptive 


thermonuclear supernova or implode to form a neutron stat25) 


The probability of collapse is highly uncertain, as it depends on the timescales of sedimenta- 
tion and crystallisation, both unconstrained at these high masses. However, the possibility of this 
novel formation channel for neutron stars is intriguing. If, upon collapse, no angular momentum 


is lost and magnetic flux is conserved, the newly born neutron star, with a magnetic field strength 
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of ~ 2 x 10!° G and a spin period of ~ 15 ms, would resemble a young pulsar. Owing to the 
gradual collapse, the neutron star would likely not receive a strong natal velocity kick. We would 
therefore expect such a neutron star population to be more concentrated to the Galactic plane. Fur- 
thermore, the proximity of ZTF J1901+1458 (41 parsecs) means it is not a rare type of object, and 
thus this formation channel would contribute appreciably to the total neutron star population. ZTF 
is currently discovering large numbers of similarly massive and rapidly variable white dwarfs. This 


enlarged sample will help us better understand the origin and fate of such objects. 
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Table 1: ZTF J1901+1458 parameters. The ephemeris Tọ, Barycentric Modified Julian Date in 


barycentric dynamical time (BMJDrgp), corresponds to a minimum in the light curve and it is 


given separately for the ZTF g-band and r-band because the latter lags behind the former by 51 


seconds. Here, Tọ corresponds to the minimum of each lightcurve. 
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Methods 


Period Detection Our search for periodicity in massive white dwarfs was part of a broader search 
for periodic variability on and around the white dwarf cooling track with ZTF, which has already 
yielded several results, including finding numerous double white dwarf binarie? B364, The targets 
were selected using Pan-STARRS (PS1) source catalogu! cross-matched with a white dwarf 
catalogu’? after imposing a photometric colour selection of (g — r) < 0.2 and (r — i) < 0.2. As 
the sensitivity of period finding depends strongly on the number of samples in the lightcurve, we 
limited the search to those targets for which 50 or more photometric 5c detections are available 
in the ZTF archival data. In order to maximise the number of epochs for each lightcurve, we 
combined data from multiple filters by computing the median magnitude in each filter, and shifting 
the g- and ?-band so that their median magnitude matched the r-band data. We used a graphics 
processing unit (GPU) implementation of the conditional entropy period finding algorithm“ We 
cross-matched our candidates with the Gaia DR2 catalogue? and visually inspected the lightcurves 
of those objects that lie below the main white dwarf cooling track in the Gaia colour-magnitude 
diagram. ZTF J1901+1458 stood out because of the high-significance detection of its short period 


and its blue and faint location in the colour-magnitude diagram (see Fig. [1). 


Magnetic Field At low magnetic field strengths, the effect of magnetic field on hydrogen transi- 
tions can be calculated as a perturbation to the zero-field electron wavefunction. The perturbation 
lifts the degeneracy in the m quantum number, and each transition is split into three Zeeman com- 
ponents, corresponding to a change in m of +1, 0 and —1. In the strong regime (above ~ 100 MG), 


the magnetic and Coulomb terms are comparable, and the wavefunction does not have a spherical 
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symmetry anymore. The perturbation method is not viable in this case, and the energies and oscilla- 
tor strengths of hydrogen transitions in this regime were calculated using variational methods? 
An important characteristic of high-field transitions is that some transitions become ‘stationary’, 
i.e. they go through a minimum or a maximum in energy and, at those field strengths, appreciable 
changes in B only yield small changes in wavelength. The magnetic broadening at these transi- 
tions is considerably reduced, and therefore the absorption features are still detectable even after 
averaging the field over the surface of the white dwarf. We considered all the bound-bound transi- 
tions tabulated in Ruder et al. 19943 that follow selection rules, and we found that the features in 
the spectrum of ZTF J19014+1458 correspond to stationary lines in a ~ 800 MG field. The wave- 
lengths of the transitions in this regime as a function of magnetic field strength are shown in the 
lower panels of Fig.|3} The identified features are numbered and the corresponding transitions are 
listed in Extended Data Table 2] where the states are labelled by their asymptotic quantum num- 
bers nlm for the zero-field states and by the numbers NM u for B — oo, where u defines the z 
parity of the states [7 = (—1)“]. A few lines show variations with phase (see Extended Data Fig [4] 
and [5). which can be attributed to the variation of the magnetic field strength across the surface. 
In particular, some of the features become narrower or broader at different times, indicating that 
the field is more or less uniform over the surface at different rotation phases. Additionally, the 
feature at ~ 4600 Å(line 7) shifts in wavelength, going as low as ~ 4500 Å, which indicates that 
the magnetic field is as low as 600 MG on some regions on the surface of the white dwarf. This is 


also confirmed by the contamination of line 10 by line 11. 
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Photometric Fitting To determine the radius and effective temperature of ZTF J1901+1458, we 
made use of the available Pan-STARRSHS photometry and the Gaia? El parallax. In addition, we 
obtained Swift UVO photometry (TOO proposal number 14380, target ID 13650). The photo- 
metric data used in the fitting is listed in Extended Data Table|[1] As the errors in the Pan-STARRS 
photometry are lower than the photometric variation, we used an error of 0.02 magnitudes instead 
to take into account the error induced by variability. In order to estimate the temperature and red- 
dening, we fitted the photometric data with the white dwarf 1-D model DA (hydrogen dominated) 
atmospheres of Tremblay et al. 201 122. The value of the radius is then estimated thanks to the paral- 
lax measurement from Gaia. In order to account for extinction, we applied reddening corrections to 
the synthetic spectra (and polynomials) using the Cardelli et al. 198953 extinction curves, available 
athttps://www.stsci.edu/. From the corrected models, we computed synthetic photom- 


etry using the pyphot package (https: //mfouesneau.github.io/docs/pyphot/). 


For the fit, we used a Levenberg-Marquardt algorithm, and the free parameters were the 
effective temperature 7;g and radius R, of the white dwarf as well as the interstellar reddening 
E(B —V ). For the reddening, we imposed a prior based on the distribution of interstellar extinction 
measured by Gaia for the stars in the same area of the sky as ZTF J1901+1458. In particular, we 
selected the stars that are within 5 degrees of the location of white dwarf and within 50 pc from 


Earth (see Extended Data Figure). 


The best fit is shown in Extended Data Fig. |1| while the left panel of Extended Data Fig. 
illustrates the corner plots. In the fit, we assumed the nominal value of the Gaia parallax for the dis- 


tance. As a change in distance would have the effect of changing the normalisation in the spectrum, 
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the error in the distance only influences the error in the stellar radius, and therefore we included the 
error in the distance in the error in the radius. The errors quoted are the statistical errors derived 
from the Monte Carlo simulation shown in Extended Data Fig.|3| They do not include possible 
modelling errors due to the effect of the magnetic field; however, as the photometric variation is 
quite small, we do not expect the magnetic field to introduce a significant error on the synthetic 


values. 


In order to verify whether our estimates are robust against the uncertainties in the spectral 
models, we employed a different set of models for hydrogen dominated atmospheres, the ones 
developed by Bohlin et al 2020“) and performed the same analysis introduced above. The results 
are shown in the right panel of Extended Data Fig. |3} the estimates of temperature, radius and 
reddening all agree within the 1c errors. Magnetic models that tackle the effects of magnetic field 
on the continuum opacities are currently unavailable, and therefore we cannot estimate the errors 


due to neglecting the field. 


Theoretical Mass-Radius Relation In order to derive the mass of the white dwarf, we computed 
the mass-radius relation employing an equation-of-state that includes Coulomb corrections to the 
pressure and energy of a degenerate Fermi gas Ed as done in Hamada and Salpeter 196 1 (H&S). 
Since ZTF J1901+1458 is so compact, general relativistic corrections are important and therefore, 
contrary to H&S, we did not integrate the Newtonian hydrostatic equilibrium equation, but rather 
the Tolman-Oppenheimer-Volkoff (TOV) equation 849) The results are shown in Fig. [4] for sev- 
eral compositions. Our calculations were performed assuming the zero-temperature limit because 


the temperature of the white dwarf is insufficient to change its structure significantly. For the 
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composition of the carbon-burning ash, we use the results from Camisassa et all For the density 
thresholds for the electron capture onto magnesium and neon) we use the values of 4 x 10? gcm-? 
and 9 x 10? g cm? respectively. For sodium?! we use 1.7 x 10? g cm-?. We estimate the tem- 
perature at a density of 1.7 x 10? g cm ? to be about 3 x 10” K; therefore, the neutrino cooling 
ZTF J1901+1458, unlike for all other known white dwarfs, is cooling dominated by the Urca 
process on ?3Ng?4P3! which makes an age determination difficult. A recent work24 simulated the 
evolution of Urca-cooling white dwarfs with the Modules for Experiments in Stellar Astrophysics 
code (MES A SE5SI showing that the Urca processes are the main cooling mechanism in the core of 


stars like ZTF J1901+1458 when they are younger than about 30 Myr. From the same models, we 


can estimate the cooling age of the white dwarf to be between 10 and 100 Myr. 


Sedimentation and collapse The solid red curve in Fig. 4] was calculated assuming that the core 
composition of the white dwarf is homogeneous and that all the sodium currently at densities 
above the threshold for electron capture has already undergone inverse beta decay, and similarly 
for magnesium and neon. As the sedimentation proceeds, much of the remaining ??Na will sink to 
densities above the threshold and undergo beta-decay, reducing the number of electrons in the star 
and reducing its radius. This means that the solid red curve in Fig. [4] will be lowered, because the 
equilibrium radius for any given mass will be smaller, and therefore the red triangle (the maximum 
mass allowed) will sink. This process can only be stopped if the core crystallises before enough 
sodium can reach the centre, as crystallisation would de facto freeze the composition gradient. 
If the star lies at the small end of the radius constraint and if at least sixty percent of the Na 


manages to sink and decay before the core crystallises, electron-capture on 24Mg would ensue and 
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the radius of the white dwarf would shrink to about 1,550 km. The central density at this point 
would be 6 x 10? g cm~, still below the threshold for electron capture on neon, but the mass of 
the white dwarf would be above the maximum mass (the red triangle in Fig. 4| would have sunk 
below the current white dwarf mass). The internal pressure would be then insufficient to support 
the star, and the star would begin to collapse, heat up and start electron capture onto neon and 


nuclear burning of oxygen. 


The possibility of collapse is highly uncertain as it depends on the timescales for crystallisa- 
tion and for the sedimentation of sodium, both currently poorly constrained, especially for such an 


extremely massive white dwarf. 


Origin of Photometric Variability The photometric variability of ZTF J1901+1458 may arise 
from either rotation or pulsations, but we believe rotational modulation is more likely. A large 
fraction of hot DQs (hot white dwarfs with carbon-dominated atmospheres) are magnetic and pho- 
tometrically variable, often with periods of several minutes? 97] Although ZTF J1901+1458 is not 
a DQ white dwarf, its other properties such as high mass, high temperature, high magnetic field, 
and rotational modulation are similar to the hot DQ white dwarfs, albeit more extreme. A useful 
case study is that of the white dwarf SDSS J142625.71+575218.3, which varies with a period of 
P ~ 418 seconds? and has a magnetic field strength of B ~ 1.2 MG! Non-radial g mode pul- 
sations were originally thought to be responsible for the photometric modulation of SDSS J1426. 
Follow-up observations revealed a tentative additional periodicity at 318 seconds™| but additional 
data ruled out the existence of this periodicity and set stringent limits on additional periodicities 


in SDSS 1426 and other variable DQ white dwarfs?8 This led Williams! to conclude that rota- 
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tional modulation is the most likely source of variability in SDSS 1426 and other hot DQ white 
dwarfs. Presently, there are no known strongly magnetic pulsating WDs, though some magnetic 


WDs could be pulsating below detectable levels. 


Furthermore, the inferred temperature and surface gravity of ZTF J1901+1458 characterise 
it as an unlikely pulsator for known instability mechanisms. The temperature of ZTF J1901+1458, 
~50,000 K, is much higher than the predicted blue edge of the ZZ-Ceti instability strip, located 
at a temperature of about 12,500-14,000 KEIEN and of the helium white dwarfs (DBV) instability 
strip, at about 30,000 (the detection of hydrogen also discourages the DBV interpretation). 
The DOV or GW Vir instability strip includes similar and higher temperatures than what we found 
for ZTF J1901+1458, but both instability mechanisms involved, the x — y mechanism for carbon 
and oxygen and the e—mechanism for hydrogen, are inefficient at such high surface gravities (see 
for example Fig. 8 of Quirion et al 201233 and Fig. 6 of Maeda & Shibahashi 201432). It is 
also possible that magnetic fields suppress gravity mode pulsations by converting gravity waves 
to Alfvén waved®"5 which is likely the cause of suppressed dipole modes in red giant stars) 
Follow-up work showed that magnetic fields greater than B ~ 0.1 MG are sufficient to suppress 
g modes of typical periods in ZZ Ceti stars) so magnetic suppression is a possibility for ZTF 


J1901+1458. 


For these reasons, we believe the variability of ZTF J1901+1458 is most likely to be caused 
by rotation rather than pulsations, though we cannot rule out either mechanism. Follow-up obser- 
vations can place more stringent limits on the presence of other non-harmonic periodicities which 


would be expected in the pulsation hypothesis. Finally, even if pulsations are the source of variabil- 
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ity in ZTF J1901+1458, it would further enhance the extraordinary nature of this star by making it 


the most massive pulsating white dwarf and the only known magnetic pulsating white dwarf. 


Data Availability Upon request, I.C. will provide the reduced photometric lightcurves and spectroscopic 
data, and available ZTF data for the object. The spectroscopic data and photometric lightcurves are also 
available in the GitHub repository https: //github.com/ilac/ZTF-J1901-1458, while ZTF data 
is accessible in the ZTF database. The astrometric and photometric data are already in the public domain, 


and they are readily accessible in the Gaia and Pan-STARSS catalogues and in the Swift database. 


Code availability We used the pyphot package (https: //mfouesneau.github.io/docs/pyphot/) 
and the corner.py package The LRIS spectra were reduced using the Lpipe pipelind™) Upon re- 


quest, I.C. will provide the code used to analyse the spectroscopic and photometric data. 
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Extended Data Figure 1: Photometric fit. The blue solid line shows the best-fitting model spec- 
trum, fitted to Pan-STARRS and Swift photometry to determine Teg, Rẹ and E(B — V). The 
synthetic photometric values (obtained from the black line) are shown in red, while the Swift val- 
ues are shown in green with 1 c error bars and the Pan-STARRS values in blue with the error that 


we chose to account for the photometric variability (0.02 mag). 
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Extended Data Table 1: Photometric data for ZTF J1901+1458. 


PS1—z PS1-2 PS1-r P$1-9g 


16.595 + 0.006 16.276 + 0.003 . 15.977 + 0.003 15.513 + 0.002 


UVOT-U UVOT-UVW1 UVOT-UVM2 UVOT-UVW2 


14.94 + 0.06 14.42 + 0.06 14.40 + 0.06 14.38 + 0.06 
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Extended Data Figure 2: Gaia E(B — V) of the closest stars. The plot shows the extinction Ag 
as measured by Gaia of the stars within 5 degrees of ZTF 19014-1458 as a function of distance and 
converted to E(B — V) assuming a reddening law with Ry — 3.1. We use the average reddening 


or the closest stars as a prior for the fitting. The error bars show 1oc errors. 
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Extended Data Table 2: Identified Balmer transitions. The numbers correspond to the identified 


transitions in Fig. land the wavelengths of the transitions at 800 MG are fron!) For 7 and 11, we 


also list the wavelength for 600 MG (left). 


line nim > n''m NMy > N'M'w A [A] 
1 | Balmer Hy 259 — 5p, 0020011 3344.23 
2 | Balmer Hy 2s => 5p, 002 0—110 3370.41 
3 | Balmer Hy 280 > 5f6 002—009 3464.11 
4 | Balmer HG 259 — 4p 002 = 007 3694.56 
5 | Balmer H8 2sọ — 4p' , 0020-16 3805.82 
6 | Balmer HB. 289 — Afj 002 005 4249.86 
7 | Balmer HB = 2s9 > Af', 002— 0—14 4517.97, 4590.9 
8 | Balmer Ha 2p. > 3d 0—10 > 0-20 5883.63 
9 | Balmer Ha 2s; — 3po 002 = 003 6615.06 
10 | Balmer Ha 2p9 > 3d 4 002 = 0—11 8525.69 
11 | Balmer Ha 2sọ — 3p_1 001 — 0—12 8537.76, 9542.62 
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Extended Data Figure 3: Corner plots. Corner plots for the photometric fitting Left: results for 
the model atmospheres of Tremblay et al. 20112 Right: results for the model atmospheres of 


Bohlin et al 202048 
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Extended Data Figure 4: Phase-resolved spectra, blue side. The LRIS phase-resolved spectra of 
ZTF J1901+1458 in the blue side. Some small variations can be observed in the spectral features 


with phase, especially in features at ~ 4600 A and at ~ 3800 A. 
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Extended Data Figure 5: Phase-resolved spectra, red side. The LRIS phase-resolved spectra of 
ZTF J1901+1458 in the red side. Some small variations can be observed in the spectral features 


with phase, in particular the feature at ~ 6620 A becomes broader and narrower with phase. 
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